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Interactions in one-dimensional (ID) electron systems are expected to cause a dynamical separa- 
tion of electronic spin and charge degrees of freedom. A promising system for experimental observa- 
tion of this non- Fermi-liquid effect consists of two quantum wires coupled via tunneling through an 
extended uniform barrier. Here we consider the minimal model of an interacting ID electron system 
exhibiting spin-charge separation and calculate the differential tunneling conductance as well as the 
density-density response function. Both quantities exhibit distinct strong features arising from spin- 
charge separation. Our analysis of these features within the minimal model neglects interactions 
between electrons of opposite chirality and applies therefore directly to chiral ID electron systems 
realized, e.g., at the edge of integer quantum-Hall systems. Physical insight gained from our results 
is useful for interpreting current experiment in quantum wires as our main conclusions still apply 
with nonchiral interactions present. In particular, we discuss the effect of charging due to applied 
voltages, and the possibility to observe spin-charge separation in a time-resolved experiment. 



PACS number(s): 73.63.Nm, 73.40.Gk, 71.10.Pm 
I. INTRODUCTION 

One-dimensional (ID) electron systems are one of the 
theoretically best-studied examples where interactions 
lead to strong correlations such that low-energy excita- 
tions cannot be described using Landau's Fermi-liquid 
concept .El As soon as electron-electron interactions are 
switched on, electron-like quasiparticles cease to exist at 
low energies, and the elementary excitations are phonon- 
like charge and spin-density fluctuations, and topological 
zero, modes. Physical quantities of such a Luttinger liq- 
uida are determined by the velocities^ v p and v a , of the 
charge and spin-density phonons, as well as additional 
velocity parameters characterizing the energy of topolog- 
ical modes. The most striking non-Fermi-liquid feature is 
exhibited by the single-electron spectral function A(k, e) 
which is basically a measure of the integrity of an electron 
as an elementary excitation in a many-body systems In 
the absence of electron-electron interactions, the spectral 
function is given by A^(k, e) = 2Tr8(e — £&) where is 
the electronic band dispersion. For an interacting sys- 
tem, the spectral function is generally broadened. How- 
ever, in a Fermi liquid, A(k 1 e) still exhibits a distinct 
single-electron-like peak, making it possible to represent 
the system of interacting electrons as a system of non- 
interacting quasiparticles that carry the same quantum 
numbers as free electrons. Such a quasiparticle peak is 
absent in the spectral function of a Luttinger liquid. Lru 
stead, a characteristic double-peak structure appearslm 
The existence of the two peaks whose energy dispersions 
follow those of the elementary charge and spin-density ex- 
citations can be interpreted as the dynamical break-up 
of the electron into two independent entities representing 
its spin and charge degrees of freedomE 

Experimental verification of spin-charge separation es- 
sentially requires>-a direct measurement of A(k,e), e.g., 
by photoemissiom or tunnelingHii3 spectroscopy. Recent 



progress in fabrication techniques has made it possible 
to create a system of two parallel quantum wires that 
are separated by a long and clean tunnel barrier Ej Uni- 
formity of tunneling between the two quantum wires, 
labeled U(pper) and L(ower), respectively, implies that 
canonical momentum is approximately conserved in a sin- 
gle tunneling event. The possibility to tune canonical 
versus kinetic momentum by an external magnetic field 
makes it possible to perform momentum-resolved tun- 
neling studies.H3 For example, in a ID Fermi liquid, res- 
onances appear in the magnetic-field dependence of the 
linear tunneling conductance whenever Fermi points from 
different wires overlapJI3E.il i.e., when for any a, a' = ± 
the parameter 
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vanishes. [Here, «u(L) denotes the ID electron density 
in the U (L) wire when no voltage is applied, B is the 
magnetic field applied perpendicular to the plane defined 
by the two wires, and d their separation.] Close to these 
resonance condition of Eq. ([!]), the differential tunneling 
conductance (DTC) as a function of voltage and mag- 
netic field can be expected to show features arising from 
spin-charge separation. This motivates the first part of 
our study, detailed in Sec. n], where we calculate the DTC 
to lowest order in perturbation theory close to the res- 
onance point corresponding to a = a' = + for a model 
ID electron system where interactions between electrons 
with opposite chirality are neglected. It turns out that 
chargingcJ in response to the applied voltage significantly 
alters tunneling characteristics. We expect this finding 
to hokL also in the realistic case of non-chiral quantum 
wires.E3 Detailed expressions are given for the location of 
four characteristic maxima in the DTC that are a mani- 
festation of spin-charge separationO 

A perturbative treatment of tunneling is valid only 
for calculating physical properties above an, in gen- 
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eral, interaction-dependentli3 energy scale. In the ab- 
sence of interactions, this scale is given by the tunneling 
strength \t\, and the regime where perturbation theory 
fails is characterized by spatial andptemporal oscillations 
in electronic correlation functions JI3 These oscillations 
result frpin coher ent electron motion between the two 
systems 113 In Sec. Ill, we provide a theoretical framework 
to treat the nonperturbajtiise-. regime with interactions, 
present using bosonizationEJO and refermionip-tionErcJ 
techniques. In analogy with previous resultsfij a char- 
acteristic length scale Li nt = irh\v p - — u CT |/|2t| emerges 
from our calculation, measuring the relative strengths 
of tunneling and interactions. On length scales shorter 
than L int , tunneling is irrelevant, i.e., it does not affect 
the electronic structure of the two wires .Ell In the op- 
posite limit of large length scales, however, we identify 
spatial oscillations in the density response with a wave 
length L t — Trhyfv p -V a /\t\ that is renormalized from its 
value irhvp/ltl in the noninteracting limit. In addition, 
we show a peculiar mode splitting to occur in the denajfaz. 
response that is similar to the one found previouslyEjEil 
for the single-electron Greens function. Besides charac- 
teristic charge- mode velocities v p ±, an additional veloc- 
ity v = 2u p _u -/(vp_ + v a ) appears in the density re- 
sponse that should be observable, in principle, in a time- 
resolved experiment. Its existence is a manifestation of 
spin-charge separation in the tunnel- coupled double-wire 
system. The naive expectation that the density response 
is sensitive just to the charge mode is satisfied only in the 
previously mentioned limit where interactions dominate 
tunneling. 

In our study of spin-charge separation in tunnel- 
coupled quantum wires, we consider a particular model 
for an interacting ID electron system. To be specific, 
the wires are assumed to be parallel to the x direction, 
located at y = and z — zu,l> respectively. The Hamil- 
tonian is given by 

H = Hq + H tu n + Hint , (2a) 
H = Y^l ahv FP ( k ~ ak Fp) + V P eV p\ c L Q/ 3 C fc<xa/3 > ( 2 t>) 
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Indices f3, f3' — U,L distinguish between the upper and 
lower wire, and a = ± between right-movers and left- 
movers. Spin quantum numbers are denoted by a, a'. 
Each wire's Fermi wave vector k-pp = ^np + ^ zp con- 
tains the effect of a magnetic fieldll B = By. A volt- 
age Vu(L) is applied to the U (L) wire, and (in general, 



aeters vp measure the resulting 
3 The wires are coupled by a 



voltage-dependent) par; 
shift of electron bands 
tunneling matrix element t chosen to be real. Fourier 
transforms of the density of spin-er electrons from the a 
branch in the [3 wire arc denoted by g qaa p- Interactions 
included in H mt are chiral, i.e., only electrons from the 
same branch (right-moving or left-moving) within each 
wire and between the two wires interact. This model 
applies directly to interacting edge channels in quantum- 
Hall bilayersLJ when each layer is at filling factor 2 and 
Zeeman splitting is negligible. Interactions between left- 
moving and right-moving electrons which are present in 
real quantum wires are not accounted for in our model. 
It turns out, however, that strong features arising from 
spin-charge separation are accurately described already 
by the chiral modeLBula For example, the redistribution 
of spectral weight due to non-chiral interacti ons .leads 
only to small additional structure in the DTC.EaO Pos- 
sibilities to go beyond the chiral model are discussed in 
Sec.lHI. 



II. PERTURB ATIVE TREATMENT 

Results presented in this section are obtained within 
lowest-order perturbation theory in the tunneling Hamil- 
tonian displayed in Eq. (^cj) . We focus on magnetic fields 
and voltages close to the resonance point where only 
right-movers tunnel. [See Eq. ([!]) for a = a' = +, and 
the inset of Fig. |l|.] In the following two subsections, we 
provide details of the calculations and give results for the 
differential tunneling conductance (DTC), respectively. 



A. Formalism 

A standard calculation! to lowest order of perturbation 
theory in H tU n yields the expression 
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for the tunneling current. Here we denoted the barrier 
length by L and the electrochemical-potential difference 
across the barrier by V = Vtj — Vl- The Matsubara 
Greens function 

J(x,t) = 

[T^U^x, t)^+l(x, ^ +l (0, 0)W+u(0, 0)) (4) 



can be calculated straightforwardlyEj using bosoniza- 
tion methods. Here we used the notation ip a +i3( x , T ) — 
J2 k e lkx Ck a+ p(T). At zero temperature, we find 



J (x, t)| t=0 
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Spin-charge separation is manifested by the occurrence 
of four algebraic singularities in the Greens function 
J{x, t). v a p = vf/3 are the velocities of spin-density exci- 
tations in the two wires, which are unaffected by interac- 
tions. The charge-density eigenmodes in the double-wire 
system have velocities 

Vpij + w p l _l |EAjl| 



v P ± = 



± 



(6) 
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that differ from the charge- mode velocities v p p = v-pp + 
Upp/iirh) of the two respective wires due to inter- wire 
interactions. Here r = tt\v p u — v p l\/\2U\jl\- Another 
consequence of inter-wire interactions is the finite expo- 
nent 9 r = l/(2vl + r 2 ). Note that, in the limit of strong 
interactions where r — > 0, the singularity in 3{x, r) asso- 
ciated with velocity disappears, and the singularity 
for w p+ changes into a pole. n++ is the resonance pa- 
rameter, defined in Eq. ([!]), that measures the distance 
of the Fermi points for right-movers in the two wires at 
zero applied voltage. 

In real quasi-lD systems, voltage-induced shifts of elec- 
tron dispersion curves, denoted here by ^u(L) e ^u(L)i de- 
pend strongly on sample details. It is for that reason 



that we treat ^u(L) as free parameters when calculating 
the differential tunneling conductance. Our general dis- 
cussion of charging effects is intended to serve as a useful 
guide to interpret experimental data. At the same time, 
we would like to point out, however, that it is possible to 
derive explicit expressions for the parameters ^u(L) for 



the model specified by Eqs. 
eralization of previous studiesE 



Application and gen- 
of charging in Luttinger 



liquids to our double-wire model system yields 

1 v pL v FV {v p+ ~v p -)v FV Vl 
Vp+v P - 2V1 + r 2 v p+ v p - V\s 

and an analogous expression for vl- Note that charging 
is strongly affected by inter-wire interactions. In particu- 
lar, in the limit of strong interactions where no charging 
would occur for separated wires, inter-wire interactions 
can drive the double-wire system into a regime where 
charging is restored. 

In the absence of inter- wire interactions, the expres- 
sion for the tunneling current given in Eq. (0) reduces to 
the familiar form 
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The functions A/3(q,s) are known exactly;Q they are the spectral functions of chiral interacting ID electron systems, 
containing right-movers only, that are parameterized by the appropriate pair of spin and charge velocities v a p , v p p : 
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Results shown in Figs. ^]-||were calculated using Eq. (||) 
for the case V\j — —Vl — V/2. To simplify the nu- 
merical calculation of the tunneling current, we have ap- 
plied Eq. (||) also to the case with inter-wire interactions 
present, which is depicted in the inset of Fig. [|. This ap- 
proximation simply amounts to neglecting the correction 
9 r to exponents in the Greens function 3(x, r) shown in 
Eq. (Q). The location of the four maxima in the DTC 
that are exhibited by the chiral modelJHI determined by 
the singularities of J~(x,t), will still be reproduced ade- 
quately as long as 8 r < 1/2. The size of spectral weight 
inbetween maxima of the DTC, however, will not be given 
correctly by this approximation. 



B. Results 

Strong features arise in the DTC from the algebraic 
singularities of the Greens function shown in Eq. (||). 
Spin-charge separation is manifested by maxima that 
form four characteristic lines as a function of magnetic 
field and voltage. This is seen, e.g., in Figf-El where we 
show the result for the band-shifting limittiTwhere ap- 



plied voltages are assumed to shift electron bands with- 
out filling them: v\j = vl = 1- The slopes of bright 
maxima are given by the inverse of the charge and spin 
velocities. [The inset of Fig. [j] shows the DTC for noHs 
interacting quantum wires in the band-shifting limit ,E3 
indicating the region (inside the square box) that is en- 
larged in the main panel where the effect of spin-charge 
separation can be observed.] In the more general case, 
however, when the quantum wires are also charged by 
applied voltages, the slope of these lines is changed. We 
have analyzed the expression for the tunneling current 
for the case of symmetric bias (Vjj = — Vl = V/2). For 
an analytic determination of characteristic equations for 
the lines of maximal DTC as a function of magnetic field 
and voltage, we use the simplistic replacement 



Ap(q,e) 
obtaining 



S(e- v a pq) + S(e - v p pq) 



K+± _ 1 + _ l-vy 

eV 2v^l 2v a \j 
k ++ _ _L I - vl _ f - vy 

eV v p l 2v a L 2u ctU 
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and corresponding expressions where U and L are 
exchanged.E3 We see that charging reduces the character- 
istic slopes of maxima in the DTC, reaching their small- 
est value in the band-filling limit = fi, = 0), which 
is shown in Fig ^|. As becomes clear from specializing 
Eqs. ( fiol) to the band-filling limit, there is at least one 
negative slope for any quadruple of velocities. This is 
very different from the band-shifting case where the tun- 
neling current vanishes in the region of-negative k++ and 
positive voltage for kinematic reasons.Blij 

For realistic quantum wires, an intermediate regime 
< ^u,^l < 1 will be realized where applied voltages 
both shift and fill electron bands. Examples are shown 
in FigJU In the main panel, the DTC for z/jj = v-^ = 0.6 is 
displayed, illustrating the fact that knowledge of charg- 
ing properties is crucial for extracting the charge and 
spin velocities from experimental data. The inset shows 
the result to be expected for finite inter-wire interactions 



(implying v p+ /v. 



'o-u 



2.56 and v p -/v a \j = 1.44), using 



fu(L) calculated from Eq. 



III. BEYOND THE PERTURB ATIVE REGIME 

Bosonization and refermionizationS are powerful 
methods enabling exact calculation of electronic corre- 
lation functions for interacting ID systems. Here we 
apply these to the tunnel-coupled quantum-wire system 
described by ihcjaodel Hamiltonian (||) , extending pre- 
vious studiesEjo In particular, we give explicit expres- 
sions for the density response function, which exhibits 
features due to spin-charge separation in the limit where 
tunneling is relevant. We note that naive straighforward 
calculation of the tujincling current within this model 
yields a zero result J2j as perfect translational invariance 
implies coherent motion of electrons between the wires 
and, hence, vanishing current flow. Experimental detec- 
tion of the tunneling current requires leads to be attached 
to the system which brea ks , tr anslational invariance and 
results in a finite current Perturbation theory actu- 
ally simulates this situation by excluding the possibility 
for electrons to tunnel twice, which is adequate onbz_.if 
the tunneling barrier is shorter than 7rhy/v-F\jv-Fi,/\t\£3 



A. Reexpressing the Hamiltonian in new variables 

The Hamiltonian H for the interacting double- wire sys- 
tem given in Eqs. (||) contains eight flavors of electrons, 
distinguished by spin, wire index, and chirality. Follow- 
ing the steps outlined in Appendix [a|, it is possible to 
rewrite H in terms of new degrees of freedom whose dy- 
namics is simpler than that of the original electrons: 



H = 



(a) 
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dx 



H 



(a) 



- dx xt ] (-idx) X 



(a) 



=2t dx i X <r> x p 



5 I dx 1X2 X3 



dx 



A 



(lla) 
(lib) 
(11c) 
(lid) 
(lie) 

(llf) 



Here, the chiral boson fields represent fluctuations in 
the total charge and spin density in the double-wire sys- 
tem; they are defined by d x $c fir = g^ a v + Q-\ah + giaV + 
Ql aL and d x ^i a '/ir = £> TqU + Q-\ a L - QiaU - Qiah, respec- 
tively. The fields Xj = Xj"^ are Majorana fermions. 
Their relation with the original electronic degrees of free- 
dom is highly nonlinear. (See Appendix |A|.) Physical ob- 
servables can be written in terms of the bosonic normal 
modes as well as the Majorana fermions. For example, 
the density of current flowing from the upper wire to the 

lower one is given by j — e/Ti ■ t ^ Q i Xi*^X3*\ an d the 
density of electrons with chirality a in the U (L) wire is 



PqU(L) — ^ PctqU(L) — 



2tt 



/ 1 \ ■ la) (1 
(+)IX2 X 3 



(«),» 



(12) 



The advantage of representing H in terms of the 
bosonic fields and the Majorana fermions is that itu 
teractions are absorbed into their respective velocities. EJ 
We find v c = (v FV + v FL + [Uvv + U Lh + 2U vh ]/irh)/2, 

V2 = V3 = («FU + V FL + [C^UU + C^LL - 2J7ul]/7T?i)/2, 

and v s = vi — V4 — (dfu + i>fl)/2. The parameter 
6 = e(z/LVi - iajVu) + ("fu + vfl)k++/2 is given in 
terms of the applied voltages and magnetic field, and 
-ffcorr contains the effect of the wires being not iden- 
tical: A c = (vpu - vfl + [Utjv ~ Cll]/tt^)/2, and 
A s = (vpu — «fl)/2. For comparison, it is useful to 
note the relations 



_ VpV + VpL Uvt 

2 irn 

V s = Vi =V A = 



V2 = V 3 



VpU + VpL _ UlJL 

2 ' irk 



(13a) 
(13b) 
(13c) 



7— c,s 



in terms of velocities defined in Sec. |[I A| . When inter- 
wire interactions are strong, i.e., 7t?i|A c /{7tjl| = r -C 1. 
we find v c v p+ and V2 — V3 » v p - . 
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B. Results for the density response function 



We consider the retarded real-time 
function 



density response 



V$,(x,t) = -iQ(r) ([Q a p(x,T),g a0 ,(O,O)]) (14) 

for the special case of identical wires (A c = A s = 0) but 
with intep.wire interaction present. This extends previ- 
ous workEH to the experimentally relevant situation in 
real double-wire systems. Hence, in the following, we 
have v c = v p+ , v s — v\ = V4 = v a , and v-i = v% = u p _. 

The density response function is then the sum "Dagi = 
L'c"'' + rjppiV^ of a contribution from the total-charge 

mode, T>i a \x,r) = 9(t) 5' (x — av p+ r) /2tt, where the 
prime denotes differentiation with respect to the argu- 
ment of the delta function, and a term originating from 
the Majorana fermions. Here ryuu = ^ll = — Tux = 
— ^lu = 1- Further simplification arises in the case 
5 = ('on resonance' according to the definition given 
in Ref. [13|) where spatial oscillations due to coherent mo- 
tion of electrons between the two wires can be expected 
to be largest. Details of our calculation can be found in 
Appendix Ib|. A characteristic length scale emerges, given 



by A, 



Trh\v n - 



\2t\, which measures the relative 



strength of tunneling and interactions. Note that u p _ 
differs from v a by the difference of intra and inter-wire 
interactions. It approaches the spin-mode velocity not 
only in the limit of weak interactions but also when in- 
teractions between and within the wires are equal. On 
length scales that are large compared to Li n t, i.e., when 
\x — avr\ > Lj nt , we find 



9(r) 



2n 



■ cos 



(2ttx\ 



x [vS(x — avr) — v p -5(x — av p -r)] /[{v p - — v)x] , (15) 

showing oscillations with wave length L t = 
irh^jVp-Va /\t\, and two propagation velocities v p _ and 
v = 2v p -v a I '(v p - + v a ). On length scales shorter 
than Lj nt , which can be of practical relevance for 
strong asymmetric interactions, the neutral-mode den- 
sity response depends on t> p _ only: (x, t) = 
9(r) 5' (x - av p -T) /2tt. 

Single- electron Greens functions iii|-Umnel- coupled ID 
electron systems have been shownE3H to exhibit three 
singularities as opposed to the two associated with spin 
and charge degrees-pf freedom in the single system. In 
one of the studies ,EJ the velocity of the additional mode 
is given by the analog of v. From our results above, 
we see that this new mode also appears in the density 
response and should—therefore be observable in a time- 
resolved experiment.tj 

Let us briefly comment on the effect of deviations from 
the ideal case considered here. When the wires are not 
identical, HctX effectively introduce interactions between 



bosonic normal modes and fictitious fermions. Addi- 
tional interaction terms appear when forward scatter- 
ing between left-movers and right-movers is included. In 
principle, these could be treated within a self-consistent 
mean-field approximation, yielding a Hamiltonian of the 
form H — ~^2 a fl~ rr with renormalized parameters. From 
our experience^ we expect such nonidealities to sup- 
press the amplitude of coherent charge oscillations. 



IV. CONCLUSIONS 

We have considered a model system for interacting 
tunnel-coupled quantum wires where interactions be- 
tween electrons of opposite chirality are neglected. Sig- 
natures of spin-charge separation are found to appear in 
the differential tunneling conductance. We give explicit 
expressions for how these features depend on spin and 
charge velocities in the wires as well as parameters mea- 
suring charging of the wires due to the applied voltage. 

Inclusion of nonchiral interactions which are present 
in real quantum wires will not affect these predictions, 
as the locations of these particular maxima in the DTC 
are captured correctly already in the minimal (chiral) 
model. However, the values of characteristic charge and 
spin velocities as well as the detailed distribution of spec- 
tral weight will certainlsi-.be changed by the presence 
of nonchiral interactions.E3 Most importantly, additional 
shadow maxima can appear which we expect to behave 
qualitatively similar to those present in the chiral model. 

Within a nonperturbative treatment of tunneling and 
interactions, we find that spin-charge separation is mani- 
fested in the density response function by the appearance 
of an additional mode which could be observed in a time- 
resolved experiment. 
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APPENDIX A: BOSONIZATION AND 
REFERMIONIZATION FORMALISM 

Here we provide details on how to rewrite the origi- 
nal interacting Hamiltonian, given by Eqs. (||), in terms 
of bosonic and fictitious fermionic degrees of freedom to 
obtain the equivalent but more easily tractable Hamilto- 
nian displayed in Eqs. (|ll]). To keep notation simple, we 
consider here only the right-moving part of the Hamilto- 
nian H, denoted by H +1 which contains the terms with 
a = + in Eqs. (^|). The remaining part of H describ- 
ing left-movers is treated analogously. First we apply 
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the bosonization procedure outlined in Ref. g0|. Using 
the representation of antibonding and bonding electron 
states in real space, given by ^ a ,± = (fa+L ± fa?+\j)/ V%, 



r(+) 



we find H+ = J2 7 = c , s H i 
with given by Eq. (11a), and 



h~ c + h~ s + m 



(+) i &(+) 



Ac = 



-| £ (*t i _*l ai _*_ ff , + * ff)+ + H.c 



(Ala) 



H 



(+) 




erer 

x + + H - c -) ( Ald ) 



Here the tunneling strength has been absorbed in dif- 
ferent Fermi wave vectors kp 7 — (fcpu + ^fl) /2 — 
■yt/v for the bonding and antibonding fields ^^-y- 
In addition to $ CiS that appeared previously, 
we have introduced new phase fields defined via 

a^cA = £ CT (<+* CT ,+ -<-* CT ,-) and s^s/tt = 

X^c - *i,_* CT ,-J- We have used the ab- 

breviations u = (wFU + «FL + [C^UU + ^LL-2l7uL]/27r?i)/2, 

34 = (CAju + C^ll - 2C/ UL )/4, A-y F = (w FU - VFf4/2, and 
A54 = (C/tju — l^ll)/2. A bosonization identitycJ relates 
the fermionic operators to corresponding bosonic 

phase fields fa^ = {fa + ifa + & [fa + 70s]) /2: 



(A2) 



Here, A/" is a normalization constant, while i^-y de- 
notes a Klein factored which acts as a ladder operator 
for particle species indexed by quantum numbers a, 7 
and obeys fermionic commutation rules. With the help 
of the bosonization identity, we ca n re write products of 
Fermi operators appearing in Eqs. (Al) entirely in terms 
of phase fields and Klein factors. For example, we find 



i(2*c-4t/C) 



(A3) 



This is how far the bosonization procedure was applied 
in Ref. ^0| for the special case of identical wires. 

We now proceed to refermionize terms in the Hamilto- 
nian containing exponentials of phase fields. To this end, 



the phase fields <&g.g are used to define new fermionic op- 
erators, essentially by applying the bosonization identity 
in reverse: 



1 



p. - g^fe.s gi*6,s(a;) 



(A4) 



Straightforward algebra shows that certain products of 
Klein factors appearing in Eq. (A2) exhibit the proper- 



ties required for Klein factors Fg^. In particular, two 
equivalent representations can be found: 



F. 



(i) 



F 7 



(2) 



F l'+ F l 
F l< + F l- 



Fr 



(1) 
(2) 



F u- F l- 

F L+ F l+ 



(A5a) 
(A5b) 



As it turns out, products of Klein f acto rs arising in 
bosonized expressions for terms in Eqs (Al) that are bi- 
linear or quadrilinear in fermion operators can be rewrit- 
ten as products of two Klein factors from the represen- 
tations introduced in Eqs. (A5). For example, we find 



- * — <r,+ *<r,+ ± c ± c 

= ifad x fa . 



i(2*B— it/v) 



Af 2 



, (A6a) 
(A6b) 



Here we applied the bosonization identity given in 
Eq. (68) of the first paper cited in Ref. Similarly, 

employing the identity Fg Fg — F^_Ff i+ , we obtain 



. tf-l pjl) e i(*c+*s-2t/C) 
; fa fa ■ 



(A7a) 
(A7b) 



Note that the sig n of the refermionized terms given in 
Eqs. ( A6b| ) and ( A7b ) is determined by the particu- 
lar arrangement of Klein factors in the initial bosonized 
form as w ell a s the correct representation of Klein fac- 
tors [ Eqs . (|A5| )1 for the new fermions. All other terms in 
Eqs. ( [AID that contain products of Fermi operators can 
be treated analogously. As a result, we obtain finally 



Hz = 



9j_ 

~2n 



4 5 {-ivd x + 2t) fa 

ipl(-id x )ipt + fa{-id x )fa } , 
ipl(—ivd x )i>s 
+fj: [4(-id x )4 + fa{-id x )fa] } , 

~ [fa -4) fa + 4) , 

d x $ c fa ~ Vte) (4 + fa) 
d x $ s fat + fa) fa - fa 



A, 
+A 



(A8a) 

(A8b) 
(A8c) 

(A8d) 



Using the representation of chiral Dirac fermions in terms 
of Majorana fermions according to fa — (xi + 1x2) 
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and ips = (X3 + iXi)/V2, we obtain Eqs. (11). We would 
like to stress the point that only the correct treatment 
of Klein factors enables the unambiguous determination 
of the velocities of the Majorana modes. Our approach 
differs from the one taken in Rcf. |T] where fictitious 
fermions were defined via phase fields that are linear com- 
binations of charge and spin-density phase fields of the 
individual quantum wires and not those of the bonding 
and antibonding states used here and in Ref. Efl. 



given in Eq. ([l2]) and specializing to the case S — A c 
A s = yields with rj VV = ?/ L l = -iTvh = -Vw = 1 



v t c {x,t) = ^ (Td x t> c (x,T)d x <s> c (o,o)) , 

VJ(x,t) = (Ti X2 (x,t)i X2 (0,0)) , 
Vj(x > r) = -i(Tx3(x,T) X3 (0,0)) . 



(B2a) 
(B2b) 

(B2c) 
(B2d) 



APPENDIX B: CALCULATION OF THE 
DENSITY RESPONSE FUNCTION 

We obtain the retarded real-timcEl density response 
function defined in Eq. (^) by generalizing methods de- 
scribed in Refs. ^ to our case of interest. Again, to 
avoid unnecessary repetition, we focus only on the case 
of right-movers, i.e., a = +. We start by considering the 
time-ordered Greens function 

T%p(x,T)=-i(TQ a p(x,T)Q a p(0,Q)]) ■ ( B1 ) 

Applying the representation of charge densities in terms 
of the free boson and fictitious- fermion degrees of freedom 



J 



For the special case considered here (5 = A c = A s = 0), 
the correlation functions 2? J and T)J (x, r) are trivial: 



VT(x,t) = —,[x- 



v c t + ie 



sgn(r)]" 



4tt 2 

?>I(x,t) = — [x-V 3 T + iesgn(T)]~~ 



(B3a) 
(B3b) 



Calculation of X>J is nontrivial, as \2 is coupled to X i 
through the term . It is therefore advantageous to 
work in the representation of the fictitious Dirac fermion 
field ip% whose Hamiltonian is given by Eq. (A8a), and to 
use the identity 



V, 



\x,t) = X - {(T^cOz,t)Vc(0,0)) + (t4(x,t)4(Q,0)) - (t4(x,t)iJ,~ c (0,0)) - (Tipz(x, r)?/>|(0, 0))| 



(4) 



The correlation functions for the field ipc appearing in Eq. (Q) can be calculated. @ After some algebra, we obtain 

( 27r "Lr) „r_\ x ~ v i T + «esgn(r) 



— l V\ 

47r g 4 



sgn(r) 



2^4 x 



— e iujf cos[k(w) X] 
2n 



(5) 



Here we used the notationEil X = x/viV2, f = 
sgn(r)[wA— r]+ie, and k(uj) — \J kvxVit 2 -I- .g|w 2 . Thew- 
dependence of k(u>) introduces the crossover scale L int = 
itTi{v2 — vi)/2t that was found previouslyEil within a dif- 
ferent approach. Inspecting our general result (||) for 
T>2(x,t) in the limit of length scales that are larger or 
smaller than L; nt and transforming to the retard ed cor - 



relation function yields the results quoted in Sec. Ill B 
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FIG. 1. Density plot of the absolute value of the differ- 
ential tunneling conductance as function of voltage V and 
magnetic field in the absence of charging effects. The main 
figure is the enlargement of the region indicated by a square 
in the inset. The origin coincides with the resonance point 
k = k++ = [see Eq. (Ilj)] where Fermi points for right-movers 
in the two wires overlap. Velocity parameters are chosen such 
that Vpu/vau = 1.5, ivl/^o-u = 2, v p t,/v a u = 2.5, and in- 
ter-wire interactions are neglected. The energy scale E s is 
arbitrary. Slopes of strong (bright) features are equal to the 
inverse of the velocity ratios given above. 
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FIG. 2. Density plot of the differential tunneling conduc- 
tance for the band-filling limit. Parameters are the same 
as in Fig. |l]. Maxima follow the slopes ±5 (l — ^p si ), 
+ and > ^~U 1 + respectively. Note 

that, for any choice of velocities, there is at least one maxi- 
mum following a negative slope. 







eV/E s 



FIG. 3. Density plot of the absolute value of the differential 
tunneling conductance for an intermediate situation where 
the applied voltage both shifts and fills electron bands. In 
the main figure, we have chosen U\j — — 0.6, in addition 
to parameters used in Fig. |l[ Note that the slopes of bright 
maxima are suppressed from their value in the band-shifting 
case. The inset shows the result for the same system with fi- 
nite inter-wire interactions Uui,/(nhv a is) = 0.25 present and 
band-shifting parameters v\j = 0.25, = 0.05 calculated 
according to Eq. (bj). 



